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We study quantum state estimation problems where the reference system with respect to which 
the state is measured should itself be treated quantum mechanically. In this situation, the difference 
between the system and the reference tends to fade. We investigate how the overlap between two 
pure quantum states can be optimally estimated, in several scenarios, and we re-visit homodyne 
detection. 
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, Today, entanglement is widely recognized as one of the characteristic feature, if not the essence, of quantum physics 
0,1^. Historically, the experimental paradigm of entanglement has long been Bell tests [^|^, that is, experiments 
that demonstrate quantum non- locality. A more modern example is quantum teleportation The beauty of 
teleportation is that a quantum state dissolves at some point of space-time and appears at another without ever 
' existing in-between. Indeed, neither the quantum system held by the receiver, nor the classical information, that 
T-H , make teleportation possible, contain any information about the original state. 

A fundamental problem tackled by quantum information Science is to characterize entangled states, and much has 
been learnt along this line of research [3. However, entanglement is also a feature of some measurements, that we will 
refer to as coherent measurements, characterized by self adjoint operators whose eigenvectors are entangled states. 
Interestingly, both aspects of entanglement, i.e. states exhibiting quantum correlations and coherent measurements, 
• are exploited in an essential way in quantum teleportation^^ 
[ — [ A property of coherent measurements, that we will be interested in here, is that they allow to measure relative 
. properties of a set of quantum systems without gaining information about the individual subsystems. In contrast, 
^n ' there is no non-trivial manner to measure a relative property of classical systems without actually measuring each 
, system and computing the relative property from the measurement outcomes. For example, given two classical arrows, 
' there is no way to find out the angle between the arrows without gaining information about the direction of each 
Qh! arrow (at least in principle, in practice one can always forget about classical information). What happens, in this 
4L ■ classical setting, is that the first direction is macroscopic enough that it serves as a reference axis, with respect to 
^ which the direction of the second arrow is measured. This measurement can be performed with very high precision 
since the second arrow is macroscopic enough that it can be considered classical. 

Every measurement on a quantum system can be thought of as a measurement of some property of this system 
with respect to a reference system. This reference system is usually that macroscopic that one can safely disregard its 
quantum properties. We here want to consider the case where both the reference system and the measured system are 
treated quantum mechanically. In a sense, we want to "quantise" the reference system of a quantum measurement. 
■ We will treat two classes of relative state estimation problems: the quantum analogue of the estimation of the angle 
' between two arrows, and homodyne measurements. In Section^ we will address the issue of optimally estimate the 
(modulus of the) scalar product between two qubit states, extending on the work of Ref. ,6]. We will consider (i) the 
situation where each qubit state is represented by identically prepared qubits, (ii) the situation where one state is 
represented by an orthogonally prepared qubit pair, (iii) we will discuss how the problem can be generalised to qudits. 
In Section IlIII we will discuss why homodyne measurements can be thought of as relative state measurements. We 
conclude in Section Hvl 
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II. RELATIVE STATE MEASUREMENTS 



A. Each state is represented by identically prepared qubits 



Consider the problem of estimating the angle between two directions p. The first direction is represented by N 
qubits identically prepared in some state Itpi) and the second direction is represented by M qubits identically prepared 
in some state \tp2)- We will assume M > N. Our aim will be to estimate at best the value of |('0i|i/'2)P- This problem 
can be thought of as the estimation of the state of a quantum system, e.g. the first one, relative to an axis which is 
it self quantum, the second system. 

We will construct a positive operator valued measure {P(x)}: 

P{x) e P{x) > 0; f dx P{x) = 

Jo 

where denotes the Hilbert space of a qubit, and 53 (5^**^+*^) denotes the space of (bounded) operators acting on 
j£«.jv+Af I denotes the 2 x 2 identity matrix. 

When the outcome P{x) comes out, the value x is guessed. As a figure of merit, we will take the mean variance: 



A{P{x)) = / dx #1 #2 Prob(x I Vi,^2) {x- I (^^11^2) (1) 



Prob(x I ipi,tp2) denotes the conditional probability to get the outcome x when a measurement is performed on 
qubits prepared in the states ^i, ip2- More technically, this variance can be re-expressed as: 



A(P(x))= j dxdg.dg^ {^f''+''\7:U9ir^<A92rPix)7:+ig,)®n+ig2M^^ (2) 

In this expression, gi and g2 represent SU(2) elements and dgi, dg2 represent the Haar measure over SU{2). tt denotes 
the natural representation and tt^ denotes the irreducible representation obtained by restriction of tt^^ onto the 
symmetric subspace of the space of N qubits, Ti.^. Actually tt^ is the spin-j irreducible representation, with N = 2j. 
IV'o) is some fiducial state. 
From any povm P{x), one can construct another povm: 

Q{x) = J dg ® ^«^^(g)*)P(a;)(^«^(5) ® 7r«*^(g)). (3) 

that achieves the same error variance as P{x). Clearly, [Q{x), Tr^{g) ® T^tiig)] — Vx, g. Using the SU(2) Clebsch- 
Gordan series, the latter commutation relation can be rewritten as 

M+N 

[Q(x), 7r+(5)] = yx,g. (4) 

k=M-N 

A nice property of the second argument of the commutator Q is that no representation appears more than once. 
Consequently (Shur's lemma), Q{x) has the following diagonal form 

M+N 

Q{x) = ^ qk{x)lk, (5) 

k=M-N 

where 1^ is the projector onto the irreducible subspace supporting tt'j^ . 
The condition that {Q{x)} should be a povm is then expressed as 

1 

dxqk{x)^l, \/k^M-N,...,M + N, (6) 



qj{x) >oyx. (7) 

The score A{{Q{x)}) can now be written as 

M+N „i 

A{{Q{x)})= / dx qk{x)Pkix), (8) 

k=M-N 
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where Pk{x) is a degree-2 polynomial in x: Pk{x) = I^x"^ — 2I^a; + I^. Explicit expressions for the quantities are 
given in Appendix IXI 

Clearly, the optimal povm is given by qk{x) — S{x — x"""), where -Pfe(a;™'") — mino<2:<iPfc(a;). Since the Pk are 
polynomials of degree 2, the x™™ are readily calculated and one finds that 

M+N ,1^2 
k=M-N 

Values ofA°P*(A^, AI) for some values of N, M are given in Table HJ 



N 


1 


1 


1 


2 


2 


7 


20 


1 


M 


1 


2 


300 


2 


3 


7 


20 


oo 


A°P^{N,M) X 10^ 


7.41 


6.94 


5.57 


6.25 


5.83 


3.29 


1.45 


5.56 



TABLE I; Minimal variance A{N,M). 



Let us now comment a bit on Eq.Q. This equation tells us that the best strategy is to use a measurement whose 
elements are projectors onto subspaces invariant under 7r'»^+^^. Clearly, the space of the N + M particles at hand also 
supports a representation of the permutation group of + M objects, Sym(A^ + M). Now it is an important result 
of representation theory that the algebra linearly generated by all unitaries tt'^^^^ (g) and the algebra generated 
by permutation operators on g^.^ commutant of each other. Consequently, they have common invariant 

subspaces 131] • In our case, this means that the elements of our povm project onto subspaces that are invariant under 
permutation of particles. We interpret this fact as follows. No preferred reference frame is available to estimate the 
angle between two directions, but it is known which particle belongs to the set indicating the first (resp. the second) 
direction. Therefore, it seems natural that the only kind of properties that can be measured are those related the 
permutations that can be carried on the particles. 

Our point is more easily illustrated in the case where A^ = 1, M — 1. The optimal povm is then made of two pieces, 
the singlet, a state which changes sign when a permutation is applied, and the triplet, which remains unchanged when 
the a permutation is applied. Thus, our measurement actually tests permutation properties of our system, on the basis 
of which a guess of the relative angle is made: the singlet representation of Sym(2) makes us guess that |('0i|^2)P = 1/3 
(the states are rather antiparallel), and the triplet representation makes us guess that |('0i|V'2)P =5/9 (the states are 
rather parallel). 



B. One state is represented by one qubit, the other by two orthogonally prepared qubits 



We now turn to the situation where one direction is specified by two anti-parallel qubits. Thus, let {"00 7 "01 } denote 
an orthonormal basis of 3i, the Hilbert space of one qubit. One direction is specified by an element of SU(2), gi 
say, and the other direction is specified by g2 € SU(2). We are now given the state 7r((7i)'*^|-(/'o, V'i)''''(52)|V'o): and 
we want (again) to estimate at best {ipo\TT{gi)*TT{g2)\ipo)- Again, we are looking for a povm {P{x) G 13(!K),P(x) > 
0, /o dxP{x) = The fi gure of merit has a form similar to the one we had before: 



A - / dx dgi dg2(^o,V'i,V'o|7r(gi)®2®7r(g2)P(a;)^(5i)®'*^^(52*)l^o,0i,V'o)(a:-|(V'ok(5i)*^(52)|V'o)|')'. (10) 



Again we can restrict to covariant measurement and assume that 

[P{x),7r'^^{g)]=0, Vx,g. (11) 

The details of the extremisation can be found in appendix^ The main result is that, perhaps surprisingly, we find 
essentially the same mean variance as in the case where each state is represented by two identically prepared qubits 
|l6j | Unfortunately, we don't have any intuition on why parallel and antiparallel pairs should perform as well or not 
for our problem. More generally, the differences between antiparallel qubit pairs and parallel qubit pairs in quantum 
estimation theory are still poorly understood ^xM- 
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C. Generalisation to qudits 

The foregoing analysis can be straightforwardly extended to qudit systems of arbitrary finite dimension d. In 
appendix we have computed the (generalisation of the) formula © in the case that N ^ M =1. We have found 

A = ^ (^ + 3)^ (12) 

We see that the mean variance decreases with das~ l/d^. The fact that this variance should decrease with d could be 
expected because when the dimension increases, the overlap between two randomly drawn states tends (on average) 
to 0, i.e. the states are increasingly orthogonal, and thus easier to estimate. We also note that the povm consists 
again on projectors onto subspaces invariant under permutations of particles, i.e. the overlap between two quantum 
states is estimated upon testing permutation properties. 



III. HOMODYNE DETECTION 



We now wish to describe how homodyne detection can be thought of as a relative state measurement. (See also 
Ref. I'^ ^ homodyne measurement 0, two e-m fields impinge on the two input ports of a balanced beam 

splitter. One of the input is generally referred to as "signal" , and the other as "reference" . We will assume that the 
signal and the reference have the same frequency and the same polarisation. A photodetector is placed at each output 
port of the beam splitter. This scheme aims at measuring a quadrature of the signal field from the difference of the 
two photocurrents read on the detectors. Let a, b denote the annihilation operators for the two input ports, and c, d 
the annihilation operators for the output ports. The observable that is actually measured by the homodyne setup is 

c* c — d* d — a*b + b* a, (13) 

where c = ^(a + 6),d = '^i'^^^)- 

The reference field is assumed to be in a coherent state \ipr) — |/3e*^), where /? and 6 are two known real numbers. 
It is also assumed that /3 is so large that 6 w (6) = /3e**, i.e. the reference is a classical field. Then the homodyne 
setup measures the observable f3{a*e^^ + ae^'^), and thus indeed corresponds to measuring a quadrature of the signal 
field. We can choose the quadrature we wish to measure upon tuning the phase 9. 

Let \4>s) = X)n V'n-^=f |vac), denote the state of the signal. It is a remarkable fact that the probability to get a given 

outcome, say K, is invariant under the transformation e^'^b, for arbitrary values of (p, or equivalently 

*ri *n in4> 

\A) = W) ^ IV'r(</>)) = |/3e^(^+^^ 1^.) = 5] V^„^|vac) -> |^,(0)) = ^ Vn— ^|vac). 

n ^ n ^ 

Thus, whatever convention we choose for the absolute phase of the e-m field, this convention does not affect in any 
manner the consistency of homodyne measurement. For example, we could choose the convention where the field is 
described by 



f.27r 



/ Tr\M^))iM'l^)\®\M^))iM^)\- (14) 

We will restrict the remaining of the discussion to the case where the state we want to measure is a coherent state 
that we will denote \a). Then, one can re- write the state (|14|l in number state basis as 



OO —drvl^ + l 1 /I n , I r,\2 



E 



J ' ' ' \k)c{kl (15) 



k=0 



where \k)c = ; "^— Ivac) denotes a state of k photons in the mode , 

Assume that the mean photon number \a\ is known, arg(a) is the quantity (phase) we like to measure. But 
with respect to what? To a reference In words, instead of thinking of a signal and a reference system, we 
can equivalently think of a Poisson distribution of qubits all in the state \tp) (x a\0) + /3|1). In this description, the 
difference between the reference and the signal has completely disappeared. If the mean number of photons | /? p is 
known and very large, then homodyne measurement turns to be an estimation problem for qubits on a circle of the 
Bloch sphere 
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m 



5 



IV. CONCLUSION 

In summary, we have considered relative state estimation problems, where the reference system is itself quantum. 
We emphasized how general this concept of relative state is and that it conveys an aspect of entanglement dual to 
the most studied quantum correlation between subsystems. More specifically, we have investigated the problem of 
estimating the overlap between two (pure) quantum states in various scenarii. 

In the case where each state is represented by identically prepared qubits, we have noticed a connection between 
optimal strategies and measurements testing permutation properties of the systems at hand. It would be interesting 
to investigate this connection further in other estimation problem. 

We have also seen antiparallel qubit pairs and parallel qubit pairs play are equivalent when used as a reference axis 
with respect to which a qubit is measured. It is an interesting open problem to provide a qualitative explanation for 
this fact. 

We have also revisited homodyne measurements, and discuss why it is a relative state measurement. 
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APPENDIX A: EVALUATIONS OF ll 

Let us start with I^. From Schur's lemma, we find that 

4° = -^trlUd^rXV-o^^l^lM)). (Al) 
dmr/tT 



■-A/ 
1 



dimTi;^^ Jsu{2) 

trlfc 
diuiH^i dimTi^ 



d5tr(lfe(^+(5)|Vo>(^oP'^^i(<?)* ® 1m)) (A2) 

(A3) 



Similarly, one shows that 

^ A- 1+ ^^((^'^ ® l^o)(^o|)(|V'o^^)(V'„^^| ® 1m+i)), (A4) 



It = — — ^tr((i, ® 1*^+2)). (A5) 

Straightforwardly, /° = (jv+i)(M+i) " Unfortunately, we were not able to find expressions as simple for ij and /|. 
However, a direct computation shows that 

jl _ ^ ir'W'™) |2|^((Af+l)/2,m-Ar/2+l/2)|2 /a^^x 

-'' ~ M + 2 ^ \^{N/2,N/2){M/2,m-N/2)\ (M/2,m-N /2){l/2,l/2) I ' \^^) 



7-2 _ _|_ 1^0, ™) |2|^((A/+2)/2,m-Af/2+l)|2 

~ M + 3 ^ \^(N/2.N/2)(M/2,m-N/2)\ \^ (M /2,m-N /2)(1.1) I ' \^'> 



where ^lj^mi){j2 mi) ~ Oj "^liii "^i! i2, "^2) denote Clebsch-Gordan coefficients. 
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APPENDIX B: GENERALISATION TO QUDITS 

The irreducible representations of SU(d) are labelled by d-uples of positive integers toi, . . . jTO^ satisfying mi > 
• ■ ■ > nT-d .81 • These c?-uples are called the highest weights of the representations. We can moreover always choose 
nid = 0. The Clebsch-Gordan series for tt^ ig) tt^^ now reads 



r+«)7r+w Tr(M + N-k,k,0,...,0). 

0<k<min{M,N} 



(Bl) 



Again, no representation appears more than once in this series, so that Q{x) assumes again a diagonal form. 

The relations HA4p - HA5| l still hold. But giving the analogue of Eqs. HA6p - IIA7| l involves dealing with S\J{d) Clebsch- 
Gordan coefhcients for d > 2, which is a heavy business. Therefore we didn't carry our analysis as far as for the 
qubit case. There are however some interesting situations where the expressions (|A4|) - (|A5|) can be calculated relatively 
easily, such as the case where = M = 1, which we will discuss now. First we need expressions for 1(1, 1, 0, . . . , 0), the 
projector onto the antisymmetric subspace of two qudits, 1(2, 0, 0, ... , 0), the projector onto the symmetric subspace 
of two qudits, and 1(3, 0, 0, ... , 0), the projector onto the symmetric subspace of three qudits. We have: 



Then, using the fact that dimH^ = {d - N + l)\/m{d - 1)! 0, we find 



1(1,1,0,.. 


.,0) 


1(2,0,.. 


.,0) 


1(3,0,.. 


.,0) 


Then, using the 


/"(1, 1,0,. 


..,0) 


/"(2,0,. 


..,0) 


/i(l,l,0,. 


..,0) 


/'(2,0,. 


..,0) 


/2(1,1,0,. 


..,0) 


/'(2,0,. 


..,0) 



k,l=l 
d 



k,l=l 
d 



kA,7n—l 



1 



■trl(l, 1,0 



(dim7i;)2 

Tr^trl(l,l,0, 
(dimn)'^ 

^tr((l(l,l,0 

dimH+ ^ ' ' 



1 



dim?^. 



tr((l(2,0. 



1 



tr((l(l,l,0 



1 



dimHi 



tr((l(2,0, 



.,0) 
.,0) 



d-l 
2d ' 

d+1 
2d ' 



• • ,0) |Vo)(V'o|)(|V'o)(V'o| <S> 1(2, 0, . . . ,0))) = 
, 0) ® \^Po){Mmo){M ® 1(2, 0, . . . ,0))) = 
••,0)® |Vo)(V'o|)(|V'o)(V'o|®l(3,0,...,0))) 
,0)®|V'o)(^Ao|)(|V'o)(V'o|® 1(3,0 



d-2 
2d(d+l)' 
d + 3 



2d{d+l)' 
d~l 



d^ + M^ + 2d' 
d+5 



d^ + 3d^ + 2d' 



From these identities, we find can compute the formula © and obtain the mean variance H12(l . 



(B2) 

(B3) 

(B4) 
(B5) 

(B6) 
(B7) 
(B8) 
(B9) 
(BIO) 
(Bll) 



APPENDIX C: THE ASYMPTOTIC LIMIT 



Suppose that one direction is specified by one qubit, and the other by an infinite number of identically prepared 
Qubits . We can thus suppose that this second direction, that we choose to call z, is known with arbitrary precision 
|9j . We can therefore imagine that the first step of our measurement consists in building a classical system that will 
serve as a z-axis. We are thus (again) looking for a povm {P{x)} satisfying the conditions 



< P{x) < 1, 



dxP{x) 



(CI) 
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As a figure of merit, we will consider 

JdgJ^ dx{^o\7T{grP{x)7T{g)\i;o){x-\{i;o\7r{gMo)\'f. (C2) 
For any povm {P{x)}, the povm whose elements are 
Qix)=J ^e-''^'P{x)e''^' 

achieves the same score. We can thus assume that P{x) is diagonal in the 2;-basis: 

/ so{x) \ 
V si{x) J ■ 

The mean error can again be written as 

A = [ dx{Io{x)x^ -2h{x)x + l2{x)). (C3) 
Jo 

Let us calculate /o(.x),/i(a;),/2(a;). One readily checks that 

Io{x) = tvP{x) J dg7r{g)\iJo){H^{gT, (C4) 
h{x) = tr[(P(x)®|V'o)(V'o|) / dg{7r'^'{gM^'){^f\n^\gr)], (C5) 
hix) = tr[(P(x)®|Vo®2)(Vo®2|) I dg{7r^\gMf'){i>f'\n^'{gr)]. (C6) 
Using Shur's lemma, we get 

Io{x) = ltTP{x)^^{so{x)+s,{x)), (C7) 
h{x) = i tr[(P(x) ® m{M)S2] = lisoix) + ^si(x)), (C8) 

i2{x) = \tv[{p{x) ® i^prn^rm = \{so{x) + ls^{x)). m 

We can then simply write 

A = dx[so{x){^ -lx+\) + - I + ^)], (CIO) 

from which we infer that the optimal povm is given by so{x) = S{x — 2/3), si{x) = 6(x — 1/3). In turn the optimal 
variance isA = l/18~ .0555. 

APPENDIX D: THE ANTIPARALLEL CASE 

The Clebsch-Gordan series for tt®'^ reads tt^ ® tt^ © tt^. The problem is now more complicated because the 
representation tt^ appears more than once. As a result, the povm elements do not have an a priori diagonal form 
anymore, but only a block-diagonal form (in the basis corresponding to the irreducible representations): 



900 (a;) loo goi(a;)loi ^ 
Qi.x) = I 9io(a;)lio 5ii(a;)lii 

533 (a;) 1,33 / 



(Dl) 

Explicit expressions for the operators Iqo, Iqij lio, lii) I33 will be given below. 
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There exist again functions Io{x), Ii{x), hix) such that we can write the average error as A = Jq{I2{x) — 2Ii{x)x + 
Io{x)x'^). Now, 



= j dg{il>o, 



Vi,Vo|(l^'®7r(5r)Q(a;)(l«^®7r(5))|Vo,^i,^o)-|(V'o7r(5)|V'o)ir% (D2) 



where i = 0,1,2. We compute these expressions expUcitly. In the following, {|0),|1)} will denote an orthonormal 
basis of the Hilbert space of one qubit. We start with Io{x). 

Io{x) = trQ(a;)(|01)(01| ® ^) = trg(a;)( J dg n'^\g)\01) {01\w^^ (g)* ® (D3) 

Due to Shur's lemma, there exists constants 70 and 72 such that / dgn^^ {g)\01) {01\'k'^^ (g)* — 70^0 +'1282- So 
(resp. 5*2) is the projector onto the antisymmetric (resp. symmetric) subspace of two-qubits. Their component (in a 
computational basis) are given in terms of Clebsch Gordan-coefEcients as 

(Sum over repeated indices is understood) The constants 70 and 72 are then easily calculated: 70 = 
(01|5'o|01)/tr5o = 1/2, 72 = (01 1 5*21 01)/ tr 5*2 = 1/6. Defining Q^.TW = {uvw\Q{x)\rst), one finds that 

Io{x) = iQ^s^ix^lK^ + AZ). (D4) 
Similarly, one finds that 

h{x) = ^tr(Q(a;) ® |0)(0|)(|01)(01| ® ^2) = IqIIT{x)TZ, (D5) 
and 

l2{x) = itr(Q(x) ® |00)(00|)(|01)(01| ® 53) = ^Q'ou{x)BZ, (D6) 

wnere - <^(ir)(2s)<^(3i) 

In the computational basis, the operators 1^ are explicitly given by 

yavw _ ^(00) ^(y) ^{lu){l/2v)M00)ilw) ,j^ry. 
y^OOjxyz - "-'(l2;)(ly)'-"(00)(l2)"-'(00) ^{Ij) ' V^''' 

/-I \uvw f-i{2m) ^(Ij) ^{lu){lv) ^(2k){lw) (T\Q\ 

lillJxi/z — {lx)(ly)^ (2m)(lz)^ (2k) ^ (Ij) ' ^ °> 

(-. \uvw _ f-,{2m) M^j) ^(l«)(li;)^(00)(l«)) /p,QN 
y^ODxyz - (lx)(ly)'^ {2m)(lz)^ {QO) (Ij) ' ' 

(lio)™r = [(loi)™^]*, (DIO) 

\uvw _ ^(Im) ^(3j) ^(lu){lv) ^{2k)(lw) (J^^^\ 
\^ii)xvz - <-'(la;)(l3/)^(lm)(lz)^(2fc) ^(3j) • \^^^) 

(D12) 

With all the information that wc have gathered, an explicit calculation can now be carried to find 
Io{x) = ^qm{x) + ^q\i{x) + ^qi3{x), (D13) 

h{x) = ^qm{x) + -^mix) - ^^(501(2;) + qio{x)) + ^933(2;), (D14) 

111 1 
Hx) = ^qooix) + j^Qnix) - -^^i^oi{x) + qio{x)) + -^q33{x)- (D15) 

(D16) 

As is obvious from the block-diagonal form of Q{x), the error A can be decomposed as A = Ai + A3. A3 = 
Jq dxq33{x){^ — i/i + i/i^) and Ai can be conveniently written as Ai = dxti Q{x)F{x), where 

Q{x)=( '^0"^^) 9oi(x)\ 
V qw{x) qii{x) I 
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and where 



fnnix) = X -\ — x^, fiifa;) — a; H — x'^ , fm(x) = fmfa;) = ^x. 

Joo\ J g 2 2' ■'"^ ' 18 6 6 ^ ' °^ 6V3 

A33 can be readily extremised, setting 933(2;) = S{x — 1/2), giving A33 = 1/36. 

The extremisation of Ai is less straightforward. If we restrict to povm's with a finite number of outcomes, then 
the solutions are of the form 

Q'{x)^w^6{x-x^)^{l + n\Xi+ 71^X2 +XIX3), i = l...I. (D17) 

where dx 'J2i=i Q^i^) — 1; where Xi, X2, X^ denote the three Pauli matrices. 
Then minimising Aqi amounts to minimise 

1 " 

2 ^tfooixt) + fiiix,) + n[{fia{xi) + foiix,)) + n'j{foo{xi) - /ii(a;,), (D18) 
1=1 

with the constraints Y^'^=i = 2, J^i^i ""^f = 0) = 1j 2, 3. For 1^ — 2, our numerical extremisation has yielded 
the following povm: 

QHx) = ^S{x~xi){l+Xi), (D19) 

Q^x) = ^S{x-X2){l-Xi), (D20) 

(D21) 

where xi — .644338... and X2 — .355662.... Very interestingly, the total optimal score Aoi + A33 equals Ap^''''"°'(2, 1). 
We also wondered whether increasing the number of outcomes for the part Aqi could decrease the overall score. 
Looking for povms with more outcome, we have found no improvement. We therefore believe that the povm ljD19(l is 
indeed optimal. 
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